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Abstract 

We present theoretical calculations of the Casimir force for Au thin films near the insulator- 
conductor transition that has been observed experimentally. The dielectric function of the Au 
thin films is described by the Drude-Smith model. The parameters needed to model the dielectric 
function such as the relaxation time, plasma frequency and the backscattering constant depend 
on the thickness of the film. The Casimir force decreases as the film thickness decreases until it 
reaches a minimum after which the force increases again. The minimum of the force coincides with 
the critical film thickness where a percolation conductor- insulator occurs. 
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I. INTRODUCTION 



The attractive force between two parallel neutral plates made of a perfect conductor is 
known as the Casimir force. This force is explained in terms of the radiation pressure due 



i] 



when boundaries 



to quantum vacuum fluctuations of the quantized electromagnetic field 
are present. Casimir 's original derivation [2] substracted the vacuum energy between the 
plates, as obtained from the sum of allowed vacuum modes, from the vacuum energy when 
the plates were absent. The resulting force per unit area is given by, 

" hC7[2 (1) 



240L 4 ' 

where L is the separation between the plates. 

In 1956 Lifshitz 3] generalized Casimir's results to real materials characterized by their di- 
electric function. In this theory, the dissipative effects associated with the radiation reaction 
of the elementary atomic dipoles composing the dielectric is balanced by the fluctuating vac- 
uum field in accordance with the fluctuation-dissipation theorem. If we consider two plates 
i = 1, 2 with different dielectric functions the Casimir force is given by 



he 
2^2 



roo p J 2 

/ QdQ / dk— (G s + G p ), (2) 
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where G s = (rj^V^ exp (2k L) — and G p = (r^r^, exp (2k L) — l) -1 . In these expres- 
sions, the factors r PjS are the reflectivities for p or s polarized light , Q is the wavevector 



component along the plates, q = oo/c and k = y q 2 + Q 2 . 

Experimentally the Casimir force has been measured by several groups 



111 ] for separations between the plates as low as 20nm [12], using different techniques such as 
force balances, atomic force microscopes and micro mechanical balances. The experimental 
results verify the Lifshitz theory to a high precision. These experiments have also pointed 
out the importance that Casimir forces can have in micro and nano systems. In particular 
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the role it plays in stiction has been considered by several authors 

The Lifshitz theory |3j , requires the knowledge of the dielectric function of the materials. 
An important issue is determining which is the correct dielectric function that is consistent in 
describing the optical properties of the materials and the measurements of the Casimir force. 
Although it may be thought that the problem is straight-forward, controversial results have 
been reported. The use of of the Drude model in Lifshitz theory has been argued to violate 



Nernst's heat theorem, while the plasma model presents no problem at all [20:]. This has been 
disputed by several authors . For example, spatial dispersion in the calculation of the Casimir 
force 



21 



221 ] has been proposed to solve the problem or a non-vanishing damping constant 



in the Drude model at zero temperature . Recently, an interesting proposal suggesting 
Johnson noise as a possible finite-size effect to be included to solve this controversy was 
suggested for a system made of two parallel wires 23j. Also, Ellingsen [24| has established 
the criteria for the non violation of Nernst's heat theorem working on the real axis rather 



than the imaginary axis. At this time the apparent vio 



an interesting issue under discussion in the literature [25 



ation of thermodynamics remains 



26]. 



Several experiments and theoretical studies have been done to study how the different 
dielectric functions influence the value of the Casimir force. For example, an early attempt 
was made by Iannuzzi that measured the Casimir force between hydrogen switchable mirrors 
(HSM) [27]. HSM's change their dielectric function dramatically when immersed in an 
hydrogen rich environment. However, this has been shown to have little effect on the Casimir 
force as the optical properties of the mirrors changed only in a narrow frequency range. These 
results were refined recently, describing the conditions needed to observe changes in the force 



Detween switchable mirrors 



Besides metals, semiconductors such as Si has been used 



29l . |30) and light modulation of the Casimir force has been measured 
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331 ] . In this 



case, light changes the carrier density of the semiconductor making it more or less metallic, 
thus changing the magnitude of the Casimir force. Changing the magnitude of the force 
by a suitable choice of materials can be achieved using for example silicon based aerogels 
3_4j that have the lowest index of refraction of any solid or with temperature changes of 



the dielectric function o 
dielectric to conductor. 



materials such as VO2 that undergo a martensitic transition from 
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371 ] . Recently, the possibility of using metamaterials has also 



been considered. The effect of having a negative index of refraction is to have a repulsive 



Casimir force 
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Up to now, gold is the most common metal to work with in Casimir force experiments. 
When comparing with theoretical calculations, tabulated data for the dielectric function is 



used 



401 ] . The films are usually thick enough to be able to use the assumption that the 



measured properties are close to those of bulk Au. At higher frequencies where quantum 
effects come into play, the dielectric function of Au can be described by an analytic model 
using a phenomenological approach that adds to the Drude model a Lorentz type response 
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function 



41] . Another approach based on time dependent density functional theory has 



been used to describe the intraband transitions of A u [42| . Tabulated dielectric data for Au 
depends on the conditions of the sample preparation 43j and the Drude parameters obtained 
from the extrapolation to low frequencies gives different values for the plasma frequency and 
damping parameter. This implies variations of up to 5% in the calculation of the Casimir 
force. 



Another important parameter affecting the value of t 
This was shown experimentally by Iannuzzi 44, 



re Casimir force is the film thickness 



451 ] when it was demonstrated that 



the Casimir attraction between a metallic plate and a metal coated sphere depended on the 
thickness of the coating. Most recently, in reference 46j different Au samples were prepared 
under similar conditions with thicknesses ranging from 120 nm to 400 nm. From measured 
ellipsometry data it was verified that the plasma frequency varies from 6.8 eV to 8.4 eV for 
this set of particular samples, changing the Casimir force a few percent. Also, the influence 
of slab thicknesses on the Casimir force was studied by Pirozhenko 37J] for slabs of finite 
thickness for several materials such as VO2 and doped silicon taking into account different 
carrier concentrations. 

The reduction of size can significantly change the physical parameters of a system such as 
the Debye temperature or the conductivity. The case of clusters is well known. For example, 



a metal 
cluster 



ic cluster can go from being an insulator to a conductor depending on the size of the 



4=71 ] . A similar phenomenon is observed in metallic thin films. As the thickness of 



the film approaches the mean free path, the Debye temperature and the conductivity show 
a sharp decrease in its values. This was shown experimentally by Kastle 48J with Au films 
whose thickness varied from 2 nm to 70 nm. Indeed, a conductor-insulator transition is 
observed function of film thickness in Au {^J. 

In this paper, we study the Casimir between Au thin films near the insulator-conductor 
transition. The dielectric function for the Au films is described by the Drude-Smith model 
and the parameters for this model are obtained from reported experimental data. 



II. DIELECTRIC FUNCTION OF AU THIN FILMS 

Thin films grown by thermal deposition, present anomalous behaviors in many of their 
physical properties, deviating from their bulk values with varying film thicknesses. In the 
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case of the conductivity, as the film thickness decreases it becomes less conductive until it 
reaches a point where it becomes insulating. The reason for this change can be explained 
as a percolation transition. During deposition, disordered gold islands or clusters start 
forming on the substrate confining the conduction electrons to these islands. As the depo- 
sition continues, the fraction of island covering the substrate increases making it possible 
for electrons to hop between islands and increasing the conductivity. Finally, more of these 
clusters are connected until a percolation threshold is reached J^. After this percolation 
transition, the conductivity of the film will increase with increasing thickness until the bulk 

n 

values are reached. Experimentally this has been measured by Walther [49] using terahertz 
time-domain spectroscopy. 

An analytic expression of the dielectric function that describes the Au films near this 
transition is the Smith-Drude model. This is a classical correction to the Drude model, 
introduced by Smith 5(3] to describe the conductivity and dielectric properties of disordered 
metals, liquid metals and recently the metal- insulator transition of thin Au films 49[]. This 
model assumes that the electrons in a metal are scattered with a probability that follows 
a Poisson distribution. In this model, the current as a function of time after an impulse 
electric field is applied is J(t) = J(0) exp (— £/r)/(n), where f(n) is a function that has the 
information on the probability that an electron was scattered n times in the time intervalfO, t] 
and r is the relaxation time. For the Drude-Smith model we have that 

n=l v 7 

After a scattering event, the electron will retain only part of its initial velocity. This is 
represented by the value of the coefficients c n . The Drude model assumes that after each 
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collision the electron has a velocity not related with the velocity before the collision 
thus, c n = for all n. 

For a material with a plasma frequency uo p and damping constant 7 = 1/r, the conduc 
tivity for the Drude-Smith model is 



47r(7 — iu) 



7 Cr, 



' (7 - iu) 1 

n=l v ' 



(4) 



The behavior predicted by the Drude-Smith model differs from that of the original Drude 
model at low frequencies. Let od be the conductivity predicted by the Drude model. At 
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zero frequency the DC current for the Drude-Smith model is cr(0) = ct.d(0)(1 + cl). In the 
extreme case of cl = — 1 the DC conductivity is zero. 

To calculate the Casimir force we need the dielectric function, which in the Drude-Smith 
case is given by 

i n l n c n 



e{oo) = 1 



00 (00 + ry) 



{ + 17)' 



(5) 
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Typically only the term n — 1 is enough to describe experimental results. If Ci = the 
Drude result is obtained and the case where c\ = — 1 corresponds to full backscattering of 
the carriers. Experimentally, it has been observed that mercury is described with c\ = —0.49 
and the value c = —0.94 corresponds to the quasicrystal AZe3.5Cw24.5i 7 612 [531 ] . 

In the case of the Au thin films, that we consider in this work, we use Walther's 
measurements using teraHertz spectroscopy, that show the crossing from insulator to con- 
ductor as the film thickness increases. The samples were prepared by thermal evaporation of 
Au and deposited on a Si substrate at room temperature and in vacuum. The conductivity 
measurements were best described by a Drude-Smith model. For films 20 nm thick the bulk 
values for the plasma frequency and damping parameter of Au were recovered. Decreasing 
the film thickness also reduces the plasma frequency until a critical thickness of d c = 6Anm 
when the conductor-insulator transition is observed. In table [I] we summarized the results 
of Walther 49( for some representative film thicknesses, as well as the value of the constant 
c used to describe the dielectric function. The damping parameter 7 does not change except 
close to the critical film thickness d c . 

In Figure 1 we present the plot of the Drude-Smith dielectric function e{iu) as a function 
of frequency for different Au film thicknesses. The vertical axis is normalized to the function 
e(iu)j) that is the Drude dielectric function (c = 0) for bulk Au with the parameters given 
by Walther We evaluated the dielectric function after rotating the frequency axis to 
the complex plane to — > iu, that is used in the calculation of the Casimir force using Lifshitz 
formula [3]. The three top curves in Fig. 1 show that as the thickness decreases the dielectric 
function decreases since, as seen from Table 1, the plasma frequency is also decreasing. 
However, at the critical thickness d c = 6.4 nm, the dielectric function is lower than the 
corresponding curve for the d = 4 nm curve in this frequency range. The lower curve (open 
circles) corresponds to the film thickness where the insulator-conductor transition occurs. 
For lower frequencies, we present in Figure 2 again the dielectric function for film thickness 
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of d = 6.4 nmand d = 4nm. For the latter the dielectric function decreases and goes to zero 
meaning a zero DC conductivity since c = — 1 for the thinnest of the films. In both Figure 
1 and 2 we have taken ujq = 10 16 s -1 . 

Now we consider the effect of the film thickness on the Casimir force, between a system 
made of a half space of Au and a Au thin film deposited on a Si substrate. The reduction 
factor n defined from Eq. (1) and Eq. (2) as n = is a convenient way to see the behavior 
of the Casimir force. To use Eq. (2) we calculate the reflection coefficients for a Au thin 
film of thickness d on top of a substrate given by 

7-01 _i_ „12„— 2(5 

' PS ~ 1 pS^ ,p\ 

T P' S ~ 1 i r 01 r 12 (O -2<5 W 
- 1 - ~T~ ' ps' ps° 

where r* J s are the Fresnel coefficients between material % and j where the subindex stands 
or vacuum, 1 is the Au thin film and 2 is the substrate. The optical length is defined as 



33,|52| 



5 = ^y/u 2 (e 1 (iu)-l) + (?k 2 . (7) 

In Figure 3, we plot the reduction factor as a function of separation using the parameters 
rom table 1. The dielectric function of the Si substrate is a Lorentz type os- 



of the Au film 

cillator model 29fl. As expected, the Casimir force decreases with decreasing film thickness. 
However, this trend changes after the critical film thickness of d = 6.4 nm where the perco- 
lation transition occurs. At this critical distance the plasma frequency attains a minimum 
value, as seen in table 1. After the critical thickness the force increases again although the 
film thickness is decreasing. The behavior at the distance d c where we see a sudden increase 



in the relaxation time and a decrease of the plasma frequency is consistent with t 



re singular 



behavior in the dielectric constant in metal-dielectric percolation transitions |54j. I55l|. This 
behavior is seen more clearly in Figure 3, where we show the reduction factor as a function 
of film thicknesses when keeping the separation between the plates fixed at L — 400 nm. 



Following previous works by Piroshenko and Lambrecht [29|, |37| , we further analyze the 
effect of the thin film using the optical length or phase factor 5 of the Au thin film [71 In 
particular, the difference between using the Drude, Plasma or Drude-Smith model can be 
seen in the behavior of the optical length with frequency. For example at the critical distance 
d c we calculate the optical length assuming a simple plasma model 8 P , the Drude model 3d 
and for the Drude-Smith model, and calculate the percent difference A = 100 x \3 — 3 p \/3 
and A = 100 x \8 — 3d\/3. This is shown in Figure 5, where we plotted the percent difference 



as a function of frequency. As expected, for high enough frequencies the plasma model is 
enough to describe the dielectric function of the Au thin film, if interband transitions are 
ignored. However, at low frequencies the percent difference between the phase factor when 
the Drude or plasma model are used becomes larger in particular for the plasma model. The 
optical length was calculated for the particular value ck/uo = 1 as in reference ([371]). 

III. CONCLUSIONS 

In this work we studied the effect that the film thickness has on the Casimir force for 
films near the critical thickness where a conductor-insulator transition occurs. To describe 
the dielectric function of the Au thin film we use the Drude-Smith model that describes 
the experimental data available in the literature. The parameters used in our calculations 
are best-fit parameters from experimental measurements. For thick films the force decreases 
with decreasing film thickness until a critical thickness is reached after which the Casimir 
force increases even with decreasing film thickness. The minimum value of the Casimir 
force is due to a decrease in the plasma frequency and an increase in the relaxation time 
when the metal-insulator transition occurs in the film. The use of Au thin films has the 
advantage that the growing techniques are well known, and experimentalist working on 
Casimir force measurements use gold coated spheres and substrates routinely. Depending 
on the film thickness the dielectric function deviates from the typical Drude and plasma 
models in accordance to the Drude-Smith model. As the thickness decreases the value of 
the backscattering constant c decreases changing the value of the DC conductivity. This 
could lead to experimental settings that further explore the role of DC conductivities and 
finite temperature controversies. 
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d(nm) ijj p x 10 15 s 1 1/7 (fs) c oja D 

20 13.19 18 1 

15 10.05 19 1 

10 6.28 19 1 

6.4 1.25 80 -0.7 0.3 

4 1.88 20 -1 



49]]). The last column 



TABLE I: Best fit parameters for the Drude-Smith model taken from Ref. ( 
shows the ratio of the DC conductivity for the Drude-Smith model to the expected DC Drude 
model conductivity. 
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FIG. 1: Dielectric function for thin films using the Drude-Smith model as a function of frequency. 
The parameters for each film thickness were taken from experimental data. The dielectric function 
is compared to the bulk Drude dielectric function reported in the experimental data of Walther 



49|. In all figures we take loq = 10 16 l/s. 



13 



0.018 




FIG. 2: In this figure we plot again the dielectric function for low frequencies of the two lower 
curves of Fig.l. This is for d = 6.4 nm and d = Anm. The dielectric function for the thinner film 
decreases We see a cross over since for d = Anm the backscattering constant is c = — 1 implying a 
zero DC conductivity and hence e(iuj) — > 0. 




FIG. 3: Reduction factor as a function of frequency for different Au film thicknesses. As the film 
thickness decreases the Casimir force also decreases until it reaches the critical thickness d = 6.4nm, 
after which the force increases even with decreasing film thickness. 
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FIG. 4: Reduction factor as a function of the Au film thickness. The reduction factor is calculated 
assuming a separation between the plates of L = 400nm and using the optical data for the Au 
films from Ref.( [43]). At the thickness of d c the percolation transition occurs, and the Casimir 
force attains a minimum. 
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FIG. 5: Percent difference between the optical length for the Drude-Smith model and the Drude 
model (solid line ) and the plasma model (crosses). The curves are calculated for a Au film of 
thickness d c . The low frequency behavior shows strong differences if the plasma or Drude model 
are used. In this figure = 10 16 l/s and ck/u>o = 1. 
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